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Abstract 

Background: In magnetic resonance electrical impedance tomography (MREIT), we 
reconstruct conductivity images using magnetic flux density data induced by 
externally injected currents. Since we extract magnetic flux density data from acquired 
MR phase images, the amount of measurement noise increases in regions of weak MR 
signals. Especially for local regions of MR signal void, there may occur excessive 
amounts of noise to deteriorate the quality of reconstructed conductivity images. In 
this paper, we propose a new conductivity image enhancement method as a 
postprocessing technique to improve the image quality. 

Methods: Within a magnetic flux density image, the amount of noise varies 
depending on the position-dependent MR signal intensity. Using the MR magnitude 
image which is always available in MREIT, we estimate noise levels of measured 
magnetic flux density data in local regions. Based on the noise estimates, we adjust the 
window size and weights of a spatial averaging filter, which is applied to reconstructed 
conductivity images. Without relying on a partial differential equation, the new method 
is fast and can be easily implemented. 

Results: Applying the novel conductivity image enhancement method to 
experimental data, we could improve the image quality to better distinguish local 
regions with different conductivity contrasts. From phantom experiments, the 
estimated conductivity values had 80% less variations inside regions of homogeneous 
objects. Reconstructed conductivity images from upper and lower abdominal regions 
of animals showed much less artifacts in local regions of weak MR signals. 

Conclusion: We developed the fast and simple method to enhance the conductivity 
image quality by adaptively adjusting the weights and window size of the spatial 
averaging filter using MR magnitude images. Since the new method is implemented as 
a postprocessing step, we suggest adopting it without or with other preprocessing 
methods for application studies where conductivity contrast is of primary concern. 

Keywords: Magnetic resonance electrical impedance tomography; Conductivity 
image; Noise estimation; Denoising method; Adaptively weighted spatial averaging 
filter 
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Background 

Electrical conductivity is a passive material property of a biological tissue or organ pro- 
viding diagnostic information of its physiological function and pathological state [1,2]. In 
magnetic resonance electrical impedance tomography (MREIT), we aim to visualize the 
internal conductivity distribution of the human body by injecting electrical current and 
measuring induced magnetic flux density data using an MRI scanner [3,4] . The MREIT 
technique uses the ^-component B z of the magnetic flux density B = (B x ,B y ,B z ) to 
recover conductivity and/or current density images [5-10]. 

In experimental MREIT studies of animal and human subjects [11-14], the quality of 
reconstructed conductivity images highly depends on the noise level in measured B z 
data and the injection current amplitude. For a given current amplitude, we can improve 
the image quality by reducing the noise level. If we reduce the current amplitude to 
avoid adverse effects such as electrical stimulations of nerve and muscle, the range of B z 
decreases proportionally and the measured B z data become more vulnerable to noise. In 
order not to deteriorate the conductivity image quality, we should reduce the noise level 
as well. 

There have been numerous studies to minimize the noise level in measured B z data 
by optimizing the data collection method including pulse sequences, RF coils, shim- 
ming, averaging, and so on. For example, the injected current nonlinear encoding (ICNE) 
method was introduced to reduce the noise level in B z data by extending the current 
injection time without increasing the scan time [15]. Various multi-echo and multi- 
coil techniques have been developed together with optimization methods to combine 
multiple signals from multiple coils and echoes [16]. 

Once we acquire /c-space data from a designed MREIT experiment, we should care- 
fully process the data to suppress the noise and enhance the image quality. This requires 
proper understanding of the noise characteristics. The noise standard deviation of B z is 
inversely proportional to the signal-to-noise ratio (SNR) of the MR magnitude image and 
the current injection time when the magnitude image SNR is higher than about 3 [17,18]. 
Since the magnitude image SNR varies over the image, the noise level of B z also varies. If 
the magnitude image SNR is much lower, for example, less than 3, then we should expect 
an excessive amount of noise in B z . This may happen when we scan an animal or human 
subject since there exist local regions of weak MR signal such as the lungs, gas-filled 
organs, and outer layers of bones. In those regions with very low magnitude image SNRs, 
there occur excessive amounts of noise in measured B z data. It is, therefore, important to 
properly suppress the noise effects in preprocessing or postprocessing steps. 

There are several denoising methods to process the acquired B z data [19-23]. These 
methods use either a PDE-based approach or a localized data processing approach and 
require a certain degree of manual adjustment of parameters. All of these methods are 
applied to measured B z data before conductivity image reconstructions as a preprocessing 
step. 

In this paper, we present a new method to suppress the noise effects in reconstructed 
conductivity images as a postprocessing step. We suggest incorporating a prior infor- 
mation from the MR magnitude image which is always available in MREIT. For current 
density imaging using measured magnetic flux density data, Joy et al. suppressed noise 
effects by setting some appropriate threshold in the MR magnitude image SNR [24]. A 
similar approach based on the morphological enabled dipole inversion (MEDI) has been 
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proposed in quantitative susceptibility map (QSM) to overcome the ill-posedness in its 
inverse problem [25]. 

The proposed method enhances the quality of reconstructed conductivity images by 
using an adaptively weighted spatial averaging filter. To determine the weights, we will 
define a distance function with respect to the noise standard deviation in measured B z 
data. Since the noise is inversely proportional to the MR magnitude image, we will incor- 
porate the magnitude image in the spatial filtering process. After describing the details 
of the proposed method, we will show how it performs with experimental data from a 
conductivity phantom and also animal subjects. 

Methods 

Noise in magnetic flux density data 

We inject current into an imaging object, of which timing is synchronized with a chosen 
MR pulse sequence. The usual current injection period is from the end of the exciting 
RF pulse to the beginning of the readout gradient. One may choose different current 
injection times and patterns depending on designed pulse sequences. For example, the 
ICNE-multi-echo pulse sequence extends the duration of current injection until the end 
of multiple read-out gradients. 

The externally injected current produces an internal magnetic flux density distribu- 
tion and its ^-component B z in a voxel results in an extra phase in the MR phase image. 
To remove systematic phase artifacts, we sequentially inject positive I + and negative I~ 
current pulses to obtain the following /c-space data: 



where M is the MR magnitude image, 8 is any systematic phase artifact, y = 26.75 x 
10 7 rad/T • s is the gyromagnetic ratio of the proton, and £2 is a field-of-view (FOV). Here, 
the superscript of S ± denotes a brief notation for <S + and S~. 
We extract the magnetic flux density B z by 



where r ± = Me^e^ 8 ^. 

The noise standard deviation of the measured magnetic flux density 5 Z is inversely pro- 
portional to the current injection time T c and the SNR of the MR magnitude image T as 



for T > 2.8 [17,18]. If we reduce the current amplitude I^, the range of B z decreases 
proportionally. To reduce the noise standard deviation, we have to increase the current 
injection time T c and the SNR T simultaneously. However, this is not possible mainly due 
to the T2 or T| decay of the MR signal. 

Conductivity image reconstruction 

The externally injected current / induces distributions of voltage u, current density J = 
(Jx>Jy>Jz)> an d magnetic flux density B = (B x ,B y ,B z ) inside the imaging object Q with its 
boundary dQ. The voltage satisfies the following partial differential equation: 




(1) 




(2) 
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V.(or(r)V«(r)) =0 in Q 
I — oVu-n — g on dQ 

where r = (x,y,z), n is the outward unit normal vector on dQ, and g denotes the 
Neumann boundary data subject to the injection current. From the Biot-Savart law, the 
magnetic flux density^ is related with the current density J = —crVu as 

B z(x) = — / ; ^ dr (5) 

4tt Jq 



Q |r-r'| 3 
where /zn is the magnetic permeability of the free space. 

There exists numerous image reconstruction algorithms to visualize the conductivity 
o or the current density J in the imaging object from measured B z data [5,9,16,26]. In 
this paper, we adopted the transversal /-substitution algorithm [27] since it differenti- 
ates the noisy B z data once and does not propagate the noise effects from one region to 
another. We first assume the imaging object with a homogeneous conductivity distribu- 
tion a H . Solving (4) with a H in place of a, we can compute the voltage, current density, 
and magnetic flux density, which are denoted as u H , ] H , and respectively. 

In MREIT experiments, we measure two B z data subject to two orthogonal injec- 
tion currents. Denoting them as B Z) \ and B z ^ the transversal /-substitution algorithm 
produces an image of the conductivity a as 

E V^{B z , n {v)-B^ n {v))-Vu^{v) 

a(r) = cr H (r) - ± n -^- 2 (6) 

E <Vi#(r),Vi#(r)> 

n=l 

where V 1 / := — |Q for a given scalar function/ and (•, •) denotes the scalar inner 
product. We used the computed values of two magnetic flux densities B^ x and B^ 2 cor- 
responding to two computed voltages and u^, respectively, for the homogeneous case 
of o H . Note that the noisy data of B Z) \ and B Zf 2 may deteriorate the quality of the recon- 
structed image of a. More details about the image reconstruction algorithm are described 
in [27]. 

Adaptively weighted spatial averaging 

To determine a neighborhood of a pixel, we define the following distance function: 
|Af(r)-Af(s)| 

D(r, s) := — , s e B r (rj(r)) (7) 

hit) 

where M(r) is the MR magnitude image, B r (rj(r)) is a neighborhood of the point r with a 
radius 77 (r), and h(r) is a function of the noise level in measured B z data. For each pixel s 
in the neighborhood of the pixel r, that is, s e B r (rj(r)), we define the following weighting 
factor w r (s): 

w r (s) := - e~ D ^ (8) 

b r 

where £ r := J2s e~ Di " r '^ is a normalization constant ensuring that 5Z s w r (s) = 1. We 
perform the adaptive spatial averaging of the reconstructed conductivity values as the 
following weighted sum: 



G W (X) = ^ W r (s)or(s). 



(9) 
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For the pixel at r, the distance D(r, s) measures the similarity of M(r) and M(s) in the 
surrounding region of r. Since the noise standard deviation of B z is inversely proportional 
to the magnitude image SNR, the weights are in effect adjusted by the noise level in mea- 
sured B z data. In (7), two parameters of rj(r) and h(r) regulate the extent of the spatial 
averaging filter. We now design the denominator h(r) and the radius rj(r) in (7) as 



h (r) oc sd Bz (r) and rj (r) oc sd Bz (r) . 



(10) 



We can choose the proportionality parameters of h(r) and r](r) in (10) depending on 
the quality of the reconstructed conductivity image to be filtered. 

The designed distance function D(r, s) preserves the conductivity value where the noise 
level of B z is low. To determine the characteristics of the weighting factor w r (s), we 
decompose the weighted conductivity cr w (r) to the true conductivity <r t (r) and the noise 
term: 

Mr) = "r(sMs) = W 'W + N * W) < U ) 

where <r t (s) and N r (&) denote the noiseless true conductivity and the noise term at s, 
respectively. The weighted conductivity cr w (r) and the true conductivity cr t (r) satisfy the 
following relation: 



\o w (r)-o t (r)\ = ± 



E e- D(r ' s) (^(s)-a f (r)+A/ r (s)) 

S6fi r (i)(r)) 



(12) 



where 



<£i(r) + £ 2 (t) 



E e _D(M) W»)-<T t (r)) 
seB r (i)(r)) 



and 



f 2 (r) = 



S€B r (l)(r)) 



-D(i,s) 



N t (s) 



(13) 



(14) 



For a small amount of noise, both h(r) and ?7(r) are small and, therefore, the relation (12) 
implies that both £i(r) and £2(1*) are small. This means that the method does no harm to 
the conductivity image when the noise level is low. 

For a large amount of noise, the first error term 6\ (r) becomes small if there was no 
significant variations of the conductivity values within the neighboring region .S r (^(r)) of 
the pixel at r. This is true only when the pixels with similar magnitude image values also 
have similar conductivity values. We will discuss implications of this restriction later. The 
second error term £2(1) is considerably reduced by the spatial averaging of the random 
noise. 



Imaging experiments 
Phantom imaging 

To evaluate the performance of the proposed method, we scanned a cylindrical phantom 
with four carbon-hydrogel electrodes (HUREV Co. Ltd, Korea). The phantom was filled 
with 0.4 S/m saline and included three cylindrical objects. The first (D\) was a TX-151 
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object (1.5 S/m), the second (D2) was an agar object (0.1 S/m), and the third (P3) was a 
TX-151 object (1.5 S/m) wrapped with an agar object (0.3 S/m) as shown in Figure 1(a). 
We placed the phantom inside the bore of our 3 T MRI scanner (Magnum 3, Medinus 
Co. Ltd., Korea). Using a custom-designed MREIT current source [28], we injected 10 
mA currents in the horizontally and vertically directions. The parameters of the spin- 
echo-based imaging sequence were as follows: repetition time Tr = 1200 ms, echo time 
= 15 ms, number of echoes Ne = 5, field of view (FOV) was 240 x 240 mm 2 , number 
of excitations (NEX) was 1, imaging matrix was 128 x 128, and total imaging time was 
10.24 min. Figure 1(b) shows the injected current pulses synchronized with RF pulses. We 
combined the acquired multiple echoes to optimize the SNR in measured B z data [16]. 
Figure 1(c) shows the measured B z image subject to the vertical current injection. 

Animal imaging 

The animals were laboratory beagles (2-3 years old, weighing 8-15 kg) without having 
any known disease. To prevent dribbling, we injected 0.1 mg/kg of atrophine sulfate. Ten 
minutes later, we anesthetized the dog with intramuscular injection of 0.2 ml/kg Tile- 
tamine and Zolazepam (Zoletil 50, Virbac, France). Twenty minutes later, we sacrificed 
it with an intravenous injection of 80 mg/kg KCL (Entobar, Hanrim Pharmacy, Korea). 
After clipping the hair, we attached four carbon-hydrogel electrodes (HUREV Co. Ltd., 
Korea) around the imaging area. The size of each electrode was 80 x 80 x 6 mm 3 . The 
procedure was approved by the Institutional Animal Care and Use Committee (IACUC). 
We injected currents in two mutually orthogonal directions between two pairs of elec- 
trodes facing each other. The injection current amplitude ranged from 5 to 10 mA. We 
adopted the ICNE pulse sequence. The imaging parameters were as follows: repetition 
time Tr = 1200 ms, echo time Te = 30 ms, FOV was 280 x 280 mm 2 , slice thickness 
was 4 mm, number of slice was 8, NEX was 6, imaging matrix was 128 x 128, and total 
imaging time was 60 min. Figure 2 shows two data sets we used in the paper. 

Results 

Figure 3(a) shows the MR magnitude image at the middle slice of the phantom. We seg- 
mented three regions of V\ (TX-151), T>2 (agar), and (TX-151 wrapped with agar). 
Depending on the noise level in measured B z data, we determined the denominator h in 
(7) to compute the distance function D(r, s). Figure 3(b) shows the denominator h in (7). 
Figure 3(c) shows the reconstructed conductivity image using the algorithm in (6). Due to 
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Figure 2 MR magnitude and magnetic flux density images for animal experiments, (a) and (b) are MR 

magnitude and B z images, respectively, of the upper abdominal region of the first beagle, (c) and (d) are MR 
magnitude and B z images, respectively, of the lower abdominal region of the second beagle. 



the short T2 relaxation time of the objects, the reconstructed conductivity values of the 
objects show strong noise defects. Figure 3(d) shows the filtered conductivity image using 
the proposed method. In the noisy regions of T>u i = 1, 2, 3, the values of the denomina- 
tor h were relatively large to produce smaller distance values in (7). Since the values of rj 
were also large in T>u more pixels in the neighboring region were included in the spatial 
averaging with relatively larger weights. 

We may quantify the effects of the filtering method by computing the amount of con- 
ductivity changes in a homogeneous region. We defined the following variance function 
to estimate how the reconstructed conductivity values varied in each region of interest 
(ROI): 

VvM) := l Vor l 2 ^ r ) ' for i= 1,2,3 (15) 

where \V{\ denotes the volume of the ROI T>{, Table 1 shows that the proposed method 
significantly improved the image quality. 

Figure 4(a) is the image of the denominator h in (7) for the upper abdominal region in 
Figure 2(a). The areas marked by the arrows have large values of h since their MR magni- 
tude image SNRs were low. The average magnitude intensity in the marked areas was 3.78, 
where the intensity of the entire image ranged from 0 to 30. Since the noise level of B z is 
inversely proportional to the MR magnitude intensity, the recovered conductivity values 
of the marked regions included a relatively large amount of noise. Figure 4(b) shows the 
reconstructed conductivity image without using any filtering method and we can observe 
severe noise effects in those regions. 
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Table 1 Variance of the reconstructed conductivity values in the ROIs using the variance 
function in (15) 







Vr> 2 (cr) 




Before filtering 


25.4 


15.7 


18.4 


After filtering 


3.39 


3.89 


4.81 



To compare the proposed method with a conventional denoising technique, we used 
the following reaction-diffusion iteration step [29] : 

V n+\ Xi y) = v n {Xiy) + a v • (^p-^j - P(y« - f) ^ 

v°(x,y) = f(x,y) 

where/ is the reconstructed conductivity image shown in Figure 4(b). Figure 4(c) shows 
the denoised conductivity image using the iteration method (16). We can see that the iter- 
ation method blurred the entire image. Figure 4(d) is the conductivity image after applying 
the adaptively weighted spatial averaging filter. We can see that the proposed method 
enhanced the conductivity image in the marked regions without affecting other regions. 
For the blurred image in Figure 4(c), we used the following parameters: iteration number = 
200, a = 0.1, and fidelity term = 0.01. These parameters were chosen to make the vari- 
ance of the conductivity values in the marked regions of the image in Figure 4(c) to be 
equal to that of the image in 4(d). 

Figure 5 shows similar results for the lower abdominal region shown in Figure 2(c). 
The average magnitude intensity of the marked regions was 1.97 for this case and the 
measured B z data were noisier than the case in Figure 4. Comparing two filtered conduc- 
tivity images in 5(c) and (d), we can see that the method proposed in this paper is clearly 
advantageous. 

Discussion and conclusion 

For animal or human subjects, their internal structures are heterogenous and there often 
exist local regions of weak MR signals. This makes the noise level in measured magnetic 
flux density data vary significantly for different pixels. Since we use the data to reconstruct 
conductivity images, noise effects are conveyed to the conductivity images with spatially 
varying image quality. Conventional filtering methods without considering this property, 
therefore, unnecessarily blur the entire conductivity image. 




S/m 




(b) 

Figure 4 Experiment results for the upper abdominal region, (a) Computed values of the denominator h 
in (7) for the upper abdominal region in Figure 2(a). (b) Conductivity image without using any filtering, 
(c) Filtered conductivity image by the iteration method in (1 6). (d) Conductivity image after applying the 
adaptively weighted spatial averaging filter. 
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Figure 5 Experiment results for the lower abdominal region, (a) Computed values of the denominator h 
in (7) for the lower abdominal region in Figure 2(c) and (d). (b) Conductivity image without using the 
filtering, (c) Filtered conductivity image by the iteration method in (1 6). (d) Conductivity image after 
applying the adaptively weighted spatial averaging filter. 



In this paper, we proposed a fast and simple method to enhance the conductivity 
image quality by utilizing the MR magnitude image. Noting that the noise level in mea- 
sured magnetic flux density data is inversely proportional to the pixel value of the MR 
magnitude image, we could adaptively adjust the weights and window size of the spa- 
tial averaging filter. In designing the filter, we used two parameters of r/(r) and h(r) in 
(7) to regulate the extent of the spatial averaging filter. Since we set their values to be 
proportional to the noise level in measured B z data, the proposed method adaptively 
changes the extent of the spatial averaging depending on the noise level at each pixel. It 
is important to properly choose the proportionality constants for rj(r) and h(r) in (10). 
Though we heuristically chose the constants in this paper, we plan to rigorously inves- 
tigate their effects on the filtered image and develop an automatic method to determine 
them. 

Since the spatial averaging is performed on neighboring pixels with similar MR mag- 
nitude image values, it may remove any useful conductivity contrast among those pixels. 
Using the proposed method, we can apply the spatial averaging to local regions of severe 
noise where we can not trust reconstructed conductivity values. We can, therefore, sup- 
press noise for the price of reduced image resolution within local regions with large 
amounts of noise. The proposed method negligibly influences the conductivity image in 
other regions with enough SNR. 

For the phantom image in Figure 3, we found that the spatial averaging recovered the 
correct conductivity values in the noisy regions since there was no conductivity contrast 
within each region. For the animal images in Figures 4 and 5, it was difficult to assert that 
the recovered conductivity values of the noisy regions were correct since there could have 
been some conductivity contrast among those pixels with similar MR magnitude image 
values. In application studies, we may analyze both unfiltered and filtered conductivity 
images if quantitative conductivity values of noisy regions are of primary concern. 

While the existing denoising methods in MREIT are applied to measured B z data before 
conductivity image reconstructions [19-23], the new method proposed in this paper is 
a postprocessing method, which can be applied to reconstructed conductivity images. 
Without requiring any image segmentation, one can easily implement the method with- 
out or with chosen preprocessing steps. We suggest adopting the proposed method in 
future experimental MREIT studies where conductivity contrast is of primary concern to 
extract diagnostic information. 
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